Spontaneous exciton condensation in ir-TiSe2: a BCS-like approach 
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TiSe2. Theoretical photoemission intensity maps have been generated by the spectral function 
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I. INTRODUCTION 

In the early 1960s, a new insulating phase was pre- 
dicted to possibly exist at low temperature in solids hav- 
ing small energy gaps. Jerome et alX published an ex- 
tended study of this phase developing a BCS-like theory 
of its ground state. However, at that time an experimen- 
tal realization of this phase was missing. 

The excitonic insulator phase may occur in a semi- 
metallic or semiconducting system exhibiting a small 
(negative respectively positive) gap. Indeed, for a 
low carrier density, the Coulomb interaction is weakly 
screened, allowing therefore bound states of holes and 
electrons, called excitons, to build up in the system. If 
the binding energy Eb of such pairs is larger than the gap 
Eg, the energy to create an exciton becomes negative, so 
that the ground state of the normal phase becomes un- 
stable with respect to the spontaneous formation of ex- 
citons. According to Jerome et alX, at low temperature, 
these excitons may condense into a macroscopic coherent 
state in a manner similar to Cooper pairs in conventional 
BCS superconductors. Kohn^ argued that exciton con- 
densation may lead to the formation of charge density 
waves (CDW) of purely electronic origin (neglecting any 
lattice distortion), characterized by an order parameter. 

1 r-TiSe2 is a layered transition-metal dichalcogenide 
exhibiting a commensurate (2x2x2) CDW"^ accompanied 
by a periodic lattice distortion below the transition tem- 
perature of Tc — 200 K. The origin of its CDW phase 
was controversial for a long time. Different scenarios 
have been proposed^'- , the best candidates being a band 
Jahn- Teller effect^ and the excitonic insulator phase. In 
2006, superconductivity has been discovered in TiSe2 
upon Cu intercalation, providing a renewed interest in 
this system^. Furthermore superconductivity also occurs 
for the pure compound under pressure^. Recently, angle- 
resolved photoemission spectroscopy (ARPES) data on 



1 r-TiSe2 were presented^. Theoretical photoemission 
intensity maps generated by the spectral function com- 
puted within the excitonic condensate phase model gave 
strong evidence for excitonic condensation in this mate- 
rial. To our knowledge, 1 T-TiSe2 is the only presently 
known candidate for a low temperature phase transition 
to the excitonic condensate state without the influence of 
any external parameters other than temperature. Indeed, 
as pressure is increased above 6 kbar on TmSeo.45Teo.55 
samples (allowing to control the gap size and thus the en- 
ergy necessary to create excitons), a transition to an insu- 
lating phase happens, whose origin can also be explained 
with exciton condensationi^. In this context Bronold and 
Fehske proposed an effective model for calculating the 
phase boundary of a pressure-induced excitonic insula- 
tor, in the spirit of a crossover from a Bosc-Einstcin to a 
BCS condensate^^. 

In this work, we present the theory from which we 
compute the spectral function used to describe photo- 
emission on TiSe2 and provide further support for the 
exciton phase scenario. In section [III we extend the 
model worked out by Jerome et al^ for one dimension 
to three dimensional and anisotropic band dispersions. 
The Green's functions of the different bands are derived 
here. Section IIIII first introduces the spectral function 
and its relation to photoemission. Then spectral weights 
and positions of the different bands are analyzed within 
this model. These theoretical results are compared to 
ARPES data of TiSe2. Finally, the chemical potential 
and discrepancies with density functional theory (DFT) 
are discussed before we conclude in section ITVl 



II. THE EXCITONIC CONDENSATE MODEL 

In this section, we present the model from which the 
spectral function describing photoemission on TiSe2 has 
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been computed. Jerome, Rice and Kohni have already 
treated in detail the case of a one-dimensional excitonic 
insulator. In their work, they consider a single valence 
band and a single conduction band, both isotropic. How- 
ever, for comparison with real experiment on the elec- 
tronic structure of TiSe2 an extension of the model to 
three dimensions with anisotropic band dispersions is re- 
quired. 



A. Description of the model 

The Hamiltonian of the model is composed of a one- 
electron part Hq and a Coulomb interaction part W. The 
one-electron part contains the dispersions of a single va- 
lence band ey{k) and of three conduction bands e],{k) 
(j- 1,2,3) 



a) Top view (perp. to k^) 



,{k)a^{k)a{k) + J2 ^c{k + Wi)h\{k)h{k). 



Here {k) and b\(k) are operators creating electrons with 
wave vector k in the valence band and with wave vector 
k + Wi in the conduction band labelled i, respectively. In 
the case of TiSe2 , we consider the valence band (mainly of 
Se 4p character) giving rise to a hole pocket centered at T 
and three conduction bands (mainly of Ti 3d character), 
equivalent by symmetry, giving rise to electron pockets 
centered at the different L points of the Brillouin zone 
(BZ) (see Fig. [l]for a sketch of high symmetry points in 
the BZ). The F point is separated from the L points by 
the three spanning vectors Wi = TL. The band disper- 
sions have been chosen of the form 



ei,(fc) 



2 "'X T "-1/ 



2my 



-f ti, cos 



/2^\ 
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+tc COS 
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which describe well the bands near their extrema as 
they are measured in ARPES experiment^. The unit 
vectors 6^ and ei±, pointing along the long and short 
axis of the ellipses, respectively, form a local in-plane 
basis for the electron pockets at the different L points. 



Thus, e^ii = where 



{Wi: 



Jiy, 



0) and 



Si±_ = u'i±/||wi_L|| where Wi± = (0,0,1) x Wi. The m„, 
rriL and nis are the effective masses of the valence band 
holes and of the conduction band electrons along the long 
and short axis of the electron pockets, respectively. The 
hopping parameters ty and tc represent the amplitudes 
of the dispersions perpendicular to the surface and k^A 
is the distance in reciprocal space between F and the A 
point. Parameters and e° are the band extrema of the 
bands. 

The interaction part W contains only the direct 
Coulomb interaction between electrons in the valence and 




FIG. 1: Schematic picture of the ir-TiSe2 bands considered 
in this model (near the Fermi energy Ep). (a) Top view of 
the BZ (perpendicular to kz). The Fermi surface has a hole 
pocket at F and three symmetry equivalent electron pockets 
at L, separated from F by the spanning vectors Wi, i = 1,2, 3. 
The side view of the BZ helps to situate the high symmetry 
points, (b) Schematic cut along the TL direction, showing 
the dispersions of the valence band (at F) and one conduction 
band (at L). 



the conduction bands 
W 



q,i 



where we have introduced partial electron density oper- 
ators 



Pa{q) = '^a'{k + q)a{k), pb.i{q,Wi) 



Y^bl{k + q)h{k). 



The above sums run over the first BZ. Finally the 
Coulomb potential reads 



with e being the dielectric function of the material. 

In fact for 1 r-TiSe2 there are three (4p-derived and 
Se related valence) bands at F near the Fermi energy. 
However, within a minimal model, we include only the 
highest lying valence band (as the valence band) in the 
calculations. Finally, the chemical potential is not ex- 
plicitly included in the model but it will be nonetheless 
discussed at the end of section IIIII 



B. Exciton physics 

In this paragraph, we introduce the formulation of 
the exciton physics in a similar way to Babichenko and 
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Kisele\*iS,. The aim of this discussion is to give a better in- 
sight into the concept of the excitons although it is not es- 
sential for the understanding of this paper. Furthermore 
the formula of this paragraph are not used for the results 
presented in this article. We start from the field opera- 
tor ipa{r) = '3*''' '^a(fc) which creates a hole in the va- 
lence band at position rand ■0b. = 
which creates an electron in the conduction band i at po- 
sition f. From these two entities, we construct the exciton 
creation operator 



E 

ki ,^2 



It is more appropriate to describe the exciton in terms of 
relative u and center of mass R coordinates. Due to the 
anisotropy of the electron pockets at L, we need to distin- 
guish the in-plane coordinates parallel and perpendicular 
to the long axis of the ellipses, Wi\\ . For simplicity, we ad- 
mit that Will is parallel to the a;-axis (for the other ellipses 
it is possible to generalize the following arguments using 
the ej|| and ei± unit vectors). Then we write u = r — r ' 
and Ra = m'^r'^/Ma + maTa/Ma with Mq, — m'^ + nia 
for a — x,y. In terms of the previously defined masses 
we have m'^ — my, since the hole pocket at F is isotropic, 
and rrix = and rUy = ms- Then, the exciton creation 
operator may be redefined as 



using the notation a{ka) instead of a{{kx,ky)) = a{k). 
The center of mass momentum Q = fci — ^2 and the 
relative momentum = maki^a/Ma + 'm'^k2_a/Ma arise 
naturally. 

At this point, we can expand the operator 6^ a in terms 
of the excitonic creation operator in reciprocal space 



Pa 



■nig 



a ] a\pa 



^(l>\{p,m)A\{Q,m) 



(1) 



where the coefficients appearing on the right hand side 
are the eigenf unctions of the hydrogen atom. In other 
words, the operator A^^{Q,Wi) creates an exciton hav- 
ing a center of mass momentum Q. The electron-hole 
bound state is described by the hydrogen state (j)x hav- 
ing the energy Ex = /xe''/8e^A^, e being the dielectric 
constant and 1/ fj, — 1/2)1^ being the reduced mass 
with = I/wq -I- According to Babichenko 

and KiseleviiSi (and generalizing to anisotropic conduc- 



tion bands), this hydrogen state obeys to 



Due to orthogonality of the hydrogen wave functions, the 
relation ([1]) can be inverted to 



A\{Q,w,) 



Et* / ^ ^ \ it/ ™Q ^ \ I ^ 



After this parenthesis which gave detail about the physics 
of excitons in the framework of our model, we now com- 
pute the equations of motion for annihilation operators. 



C. Equations of motion for the Green's functions 

With the help of the Hamilonian H = Hq + W, we com- 
pute the equation of motion for our electron annihilation 
operators 

d 

"^gl^iP'^) = Hp,t),H] = ey{p}a{p,t) 
+ J2 Vc{q}a{p + g, t)bl{k, t)h{k - q, t), 

q,k,i 

d 

+ Vc{q)a^ {k + q, t)a{k, t)b,{p + q, t). (2) 

q.k 

We now introduce Green's functions for the valence and 
the conduction bands 

Gy{k,t,t') = {-i){Ta{k,t)a){k,t')), 
Gl{k,t,t') = i-i){Tb,{k,t)blik,t'),) 

where we used the time ordering operator T. Their equa- 
tions of motion are derived directly from equations ^ 

E yM{Ta{p + q, t)b\{k, tMk - q, t)a^ {p, t')), 

q,k,i 

- 4(p + ^.)) GM t') = S{t - 1') 

~i J2 Veiq){Ta\k + q, t)a{k, t)h{p + g, t)b\{p, t')). 

Using Wick's theorem we simplify the four-operator aver- 
ages (...) by neglecting correlations, i.e. keeping only the 
lowest order terms. The calculation is similar for both 
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Green's functions. We get three two-operator contribu- 
tions, namely 

{Ta{p + q, t)bl{k, t)b,{k-q, t)a\p, t')) = 
{aip + q, t)h\{k, t)){Th{k - t)a){p, t')) 
^{a{p + q, t)h{k - q, t)) {Th\{k, t)a\p, t')) 
-{Taip + q,t)a\p,t')){b,{k~q,t)bl{k,t)), 

out of which only the first one remains (the second one is 
zero and the last one is a Hartree term which we consider 
as already included in the measured dispersions) so that 
the corresponding equations of motion become 

E Vciq){bUk, t)a{p + q, t)){TbS - q, t)a^ {p, t')), 

q,k,i 

(^§1 " 4(p + ^o) t, t') « 5{t ~ t') 

-I J2 {h ip + q, t)a^ [k + q, t)){Ta{k, t)bl {p, t')). 

(3) 

At this point, to go further we need to introduce the 
concept of the condensate phase. 

D. Condensate phase 

When the energy gap is smaller than the exciton bind- 
ing energy, the energy necessary to create an exciton be- 
comes negative and the normal ground state becomes 
unstable towards the spontaneous formation of excitons. 
Once temperature is low enough, these excitons may con- 
dense into a macroscopic coherent state analogous to that 
of Cooper pairs in the BCS theory of superconductivity. 

The first average on the right hand side of equations 
([3]), for k = p + q, can be expressed by the excitonic 
creation operator (equation ([1])) 

{bl{p + q, t)a{p + g, t)) = -^aIO, W^){A{{p + q, m)) 

X 

« q^*o{0,m){Al{p + q,m)). 

At sufficiently low temperature, only the lowest lying ex- 
citonic level is populated. 

By analogy with the BCS theory, we can identify the 
average in this last equation with anomalous Green's 
functions, after appropriate variable substitutions. These 
new functions are defined as follows 

F,(fc,i,i') = i~i){Tb,{k,t)a^{k,t')), 

F}{k,t,t') = (-z){Ta{k,t)b\{k,t'))- (4) 

They describe the scattering of a valence electron into 
the conduction band or inversely. Pushing further the 



analogy, we introduce the order parameter describing 
the condensate 

^,{p) = -lY,VMF}{p + q.t,t)) 
9 

= Y.^M{b\{p + <i,t)a{p + q,t)) (5) 
q 

« Y,yM^i{^.^,){Ai{p+q,w,)) 

q 

(here the anomalous Green's function definition ^ and 
equation ^ have been used). It quantifies the intensity 
of exciton formation between the valence band and the 
conduction band labelled i. Moreover, it characterizes 
the state of the system in the sense that, when the order 
parameter is different than zero, exciton condensation 
drives the system into the CDW phase (see section Hill 
for further discussion). 

E. The Green's function of the valence band 

With the help of these new elements, we can go back 
to equation In the particular case of the Green's 

function of the valence band, we obtain 

(»^-e.(p))a(p,t,0= 5{t~t') 

-J2 A,(p-f-7) FMt,t').i6) 

i 

In order to solve this equation for G„, we need to find a 
similar expression for the anomalous Green's function by 
computing its equation of motion. This procedure results 
in the following relation 

(^z^ -e^(p + «;,)) Fip,t,t') = ~A,ip)G4p,t,t') (7) 

where we have again identified the order parameter. 

Converting the time-dependence into a (imaginary) 
frequency z dependence with a Fourier transform allows 
us to solve the system of equations given by (O and (O 
for the Green's function of the valence band 

GMz)=(z-e.{p-)-J2^^0l—] . (8) 

F. The Green's function of the conduction band 

Calculating the Green's function of the conduction 
band involves a treatment similar to that of the valence 
band. From equation ([3| and with definitions jl]) and ^ 
we get 

(^^^^-elip + m)^ Gl(p,t,t') ^5{t-t') 

-A*{p)Fj{p,t,t').{9) 
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The equation of motion of is obtained with help of 
equation ^ and Wick's theorem 

(»|-ea(p))^;^(p,t,0 

= -tJ2 VM{Ta{p + q, t)b]{k, t)b,{k - g, t)b\{p, t')) 

k,q,3 

yc{q){aip + g, t)b]ip + q, t)){Tb,{p, t)b\(p, t')). 

(10) 
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The averages on the right hand side bring into play three 
b operators and present an off-diagonal term mixing bi 
with fej operators. When j — i, the last average lets ap- 
pear the Green's function Gl while i ^ j terms involve 
new Green's functions representing the scattering of an 
electron from one conduction band to another one (usu- 
ally called multivalley scattering) 

H,,{k,t,t') = {-i){Tb,{k,t)b]{k,t')). 
Their equation of motion reads 

(^*^ - el{p + m)^ H,,{k,t,t') = -A*ip)F]ip,t,t'){U) 

Thus, with the help of the definition of the order pa- 
rameter A, replacing this last definition into (fTO|) results 
in 



i^^-eM^F^{p,t,t') 



= -A,ip)Glip, t, t') - ^Am^iP, t, t')- (12) 

Equations pT|) and (fT^ together build a system of 
equations which can be solved with respect to G*, pro- 
viding us with the following expression after a Fourier 
transform to frequency space 

G^(p + Wi,z) = z - (p + Wj) 



|A,(p-)P 
-ei(p+Wj) , 



(13) 



III. RESULTS AND DISCUSSIONS 

A. The spectral function 

In the context of photoemission, the spectral function 
A{p, ri) plays a central role. It is directly proportional to 
the imaginary part of the Green's function and in the case 
of one-electron Green's functions as defined in section Ull 
it describes the one-electron removal spectrum. 



For the excitonic condensate model, we distinguish the 
spectral function of the valence band 

A^{p,n) = --lm[G4p,n + i5)] 

TT 

(d is here an infinitesimal real quantity) and that of the 
conduction band 

Al{p + w„n) = ~-lm[Gi{p,n + iS)]. 

TT 

To simplify further calculations, we rewrite the Green's 
functions (equation ^ and (US])) in the following forms 



(14) 



- \A,{me.^3^\{z^e^{p-rw^))) 

(15) 

(sym is the permutation symbol). The denominator 2?, 
common to all Green's functions, is 

V{p,z) = {z-e,{p))\{{z-e'Ap + Wi)) 



E 



= lliz-n^m 



(16) 



(here the index a refers to the four zeros of the denomina- 
tor V, while the other indices refer to the three conduc- 
tion bands). In the last line, the denominator is factor- 
ized in terms involving its four (real) zeros Qa (p) (which 
are implicitly functions of the order parameter A). These 
zeros can be calculated exactly. However their analytical 
forms are too long to be written here. 

This allows us to break apart the Green's functions 
(|14p and (|15p into rational expressions with minimal de- 
nominators, so that we can use Sokhotsky's formula 



1 



1 



= V iTT5{x — xo) 



X — + l€ X — Xq 

(V denotes the principal part) and write the spectral 
functions in rather simple forms 



A4p,n) = ^p^(p)5(f]-r!„(p)). 



A^ip + w,,n) 



4 

E 



p^^ip)SiQ~n^ip)), 



where the weights Pa{p) (which are also implicitly func- 
tions of the order parameter A) associated to the poles 
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3) Normal phase c) 



~ e^p)) UnM^- - 'ciP + ^n)) 



rand L 

A=0 meV 



|A„(p)| 



Until now, to ensure the generality of the theory, we al- 
ways kept the fc-dependence of the order parameter A(fc). 
However, in the practical analysis which will follow, we 
will use a fc-independent order parameter estimated from 
experiment. 



B. The spectral function and photoemission 

Within the sudden approximation, the contributions 
to the photoemission intensity are the spectral function, 
the matrix elements and the Fermi-Dirac distribution. In 
this paper we concentrate on the spectral function, being 
well established in the previous subsection. 

We now choose to fix the parameters that describe the 
band dispersions. They will take the valuesi^ determined 
from our previous ARPES study^. From the spectral 
functions, we can derive a first important information. 
The zeros of their denominator V are the poles of the 
Green's functions and therefore represent the renormal- 
ized electronic band positions in the system. Noting that 
a zero order parameter A in equation (|16p results in band 
positions that are not renormalized, one realizes that A 
is a good indicator of the strength of this renormaliza- 
tion. Moreover, since this denominator is the same for 
both the valence and the conduction bands, we see that 
the valence band at F is backfolded at L and that the 
conduction bands at L are backfolded at F (as expected 
from the electron-hole coupling) . This is a first indication 
of the CDW phase in the system. The situation for the 
band positions (not considering their spectral weight) is 
illustrated in Fig. O Part (a) depicts a cut through the 
Fermi surface (FS) (around the F and L points) in the 
normal phase, composed of the valence band hole pocket 
(red) and three (symmetry equivalent) electron pockets 
(blue). In the CDW phase, which is characterized by an 
order parameter A ^ meV, this FS changes into that 
of Fig. [2] (b) . Via the electron-hole coupling the L-points 
become equivalent to F (not yet considering the spectral 
weights) and all three conduction bands are backfolded 
onto the valence band. In parallel, in Fig. [2] (c), (d) and 
(e) we show the associated dispersions. To facilitate the 
comparison, in the normal phase, we superimpose the va- 
lence band (wi located at F, red continuous line) and the 
conduction bands (blue continuous line for ci and blue 
dashed line for C2,C3, normally located at the L-points 
but shifted here to F). For the CDW phase, we distin- 
guish two cases, one with a low value of the order param- 
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FIG. 2: (Color online) Schematic picture of the band posi- 
tions in the model (near Ef) in the normal phase and in the 
CDW phase. Wave vectors are expressed in a multiple of FA/, 
(a) In the normal phase (A = meV), the Fermi surface com- 
posed of the valence band at F (in red) and three symmetry 
equivalent conduction bands at L (in blue), (b) In the CDW 
phase (A 7^ meV), F becomes equivalent to L. The electron 
pockets at L, backfolded to F, produce "flower"-like Fermi 
surfaces at each newly equivalent high symmetry point, (c) 
Dispersions calculated parallel to TM (see Fig. 1(a)) around 
F and parallel to ALi around the three L-points in the normal 
phase plotted on the same graph (the minima of the different 
conduction bands ci,C2,C3 have been displaced from the L- 
points to r on the graph). (d),(e) Dispersions around F and 
along FAf in the CDW phase for A = 20 meV respectively 
A — 100. In the CDW phase, F and L become equivalent 
concerning the dispersions. 



eter (Fig. [2](d), A = 20 meV) and one with a high value 
of the order parameter (Fig. [2](e), A = 100 meV). These 
values are reasonable in comparison with experiment and 
help to understand how the CDW transition settles in. 
Once the order parameter increases to a non-zero value, 
there is a strong change in the band dispersions. The va- 
lence band vi and the conduction band C3 split, opening 
a gap between them. As the order parameter increases to 
A = 100 meV, vi and C3 repell each other further, while 
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FIG. 3: (Color online) Extrema of the renormalized bands as 
a function of the order parameter A (left). Position of these 
extrema on the band dispersions (right). 

ci and C2 stay at their original positions. 

We now turn to the discussion of the extrema of the 
bands. Fig. [3] presents the band extrema as a function 
of the order parameter A. We see that, except for small 
values of the order parameter, they display a linear be- 
haviour (this can be shown analytically from the denomi- 
nator in equation (jl6p for k exactly at F and L where the 
three original conduction bands have the same energy). 
Among the three conduction bands, only one (cs) shifts 
away from the Fermi energy. Its two minima (blue line) 
also increase linearly. Compared to the conduction band 
C3, the valence band maximum (vi) follows the inverse 
behaviour, thereby opening a gap below Ep- 

It is important to realize that when considering three 
conduction bands instead of one, the system remains in a 
semimctallic state at low temperature rather than evolv- 
ing into an insulating state, since the gap opens below 
the Fermi energy. Therefore, strictly speaking, the de- 
nomination excitonic insulator phase in this context is 
misleading and we rather adopted the expression exciton 
condensate phase. 

Besides the position of the bands the spectral function 
contains an additional crucial information, namely, the 
spectral weight (SW) carried by each band in the process 
of one-electron removal probed by photoemission. It is 
related to the numerator of the spectral function. We 
now add this feature to the previous figure and obtain 
Fig. m where the SW of the bands is indicated in gray 
scale. The evolution from the normal state (Fig. IHa)) to 
the CDW state with an order parameter of 20 meV (Fig. 
\Mh)) and 100 meV (Fig. gjc)) is shown. 

In Fig. O we focus on the SW of the bands at F and 
L (SW of the conduction band C2 is not represented here 
since it is exactly for every k along the FM and AL 
directions). Graphs [n^a) and (b) display the SW of the 
valence band (vi) at F and L, respectively, for different 
values of the order parameter A. Graphs [He), (d) and 
[5]Je) , (f ) show the SW for bands C3 and ci , respectively. 

We immediately see (Fig. S]) that, with respect to the 
SW, the backfolding is in fact incomplete even at a large 
value of the order parameter. Indeed, in the CDW phase 
with an order parameter of A = 100 meV (Fig. Hljc)), at 
F, the original valence band conserves 40% of its (normal 
phase) SW (Fig. [51^ a)), while the backfolded conduction 



band (ca) carries a SW of 60% (Fig. [^c)). At L, the sit- 
uation is more complicated, since three bands (wi, ci, C3) 
share now the SW. The original conduction band (ci) 
keeps a minimum of 67% of SW (Fig. [DJf)), while the 
other two backfolded bands, namely a symmetry equiv- 
alent conduction band (C3) and the valence band, divide 
among themselves the remaining 33% (Fig. [DJb), (d)). 

We also present the graphs for A = 20 meV. We see 
that there is a large SW loss (more than 80% a,t k — 
A^^) in the valence band, even larger than that for 
A = 100 nicV (Fig. [5ja)). What happens can be seen as 
follows. If we observe the graphs in Figs. 111(a), (b) and 
(c) close to F (i.e., the left panels) as a function of A, we 
see that at A = 20 meV the shape of the valence band 
tries to stay the same as for A = meV. This is achieved 
by a reduced SW of vi and an increased SW of C3. At 
A = 100 meV the dispersions are sufficiently different 
from the non-renormalized ones to carry more SW. In 
other words, for small values of the order parameter, the 
SW will be distributed along the parts of the dispersions 
corresponding mainly to the non-renormalized ones. A 
similar situation happens at L, as can be seen in Fig. [J] 
between bands ui, ci and C3. But this time, the original 
band (in the normal phase) is the conduction band ci, 
so that for small values of the order parameter, its SW 
is shared between vi, c\ and C3. In Fig. [H the dashed 
lines indicate a conduction band (ci) backfolded to F, 
which has a small non-zero SW as shown in Graph [5Je) 
(it is less than 1% for the values of the order parameter 
considered here). 

It should be noted that from photoemission data it is 
difficult to extract information concerning (thermally oc- 
cupied) states above the Fermi energy (set to eV here) , 
so that SW of band C3 is hardly measured in experimenli^. 



C. Comparison with experiment 

To emphasize the good agreement between our model 
and experiment, we further analyze experimental ARPES 
intensity maps^ in the light of the previous discussions 
(calculated intensity maps are not reproduced hereafter, 
see^ for more details). The data was collected at the 
Swiss Light Source with a photon energy of 31 eV on 
TiSe2 samples (at this photon energy, the normal emis- 
sion spectra correspond to states located close to the F 
point). 

Fig. [S] presents comparisons between the theoretical 
(left) and experimental (right) electronic structures at F. 
The experimental intensity maps at T = 250K (Fig. [5] 
(a)) and T — 65K (Fig. [6](b)) are compared to calculated 
bands with A = 25 meV and A = 75 meV, respectively. 
At T = 250K on the experimental side (Fig. [6l right), the 
situation is more complicated than in our model. Indeed, 
there are three Se4j)-derived valence bands, out of which 
two (black lines), are not considered in our model. The 
dashed white line corresponds to the valence band v\ of 
the model, which suffers already SW loss at T = 250K. It 
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FIG. 4: Band dispersions with their corresponding spectral 
weight at V and L, along TM and AL directions, respectively. 
Graph (a) describe the normal phase (A = meV), (b) the 
CDW phase with moderate excitonic efltects (A = 20 meV) 
and (c) the CDW phase with strong excitonic effects (A = 100 
meV). The dashed lines indicate a band (ci) having a small 
non-zero SW (see text). 



flattens at its maximum and deviates from tlie parabolic 
shape of the normal phase dispersion (this is clear from 
an energy distribution curve taken at k\^\ = A^^, not 
shown here). On the theoretical side (Fig. [5] (a), left), 
the dispersion reproduces well the experiment when con- 
sidering an order parameter of A = 25 meV. However, 
at this temperature, the system should be in the normal 
phase. Nevertheless, as in high temperature supercon- 
ductors, above the critical temperature, we expect fluc- 
tuations to persist well above Tc in the excitonic con- 
densate phasei^ii^. Thus, this non-zero order parameter 
above Tc may be understood in terms of fluctuations. 
According to the left graph of Fig. [B] (a), the bottom 
of the backfolded conduction band C3 appears just be- 
low Ep. In parallel, the experiment shows small humps 
in momentum distribution curves near Ep (not seen in 
the false colour map here). These can be attributed to 



FIG. 5: (Color online) Spectral weights of the bands at F and 
L, along FM and AL directions, respectively, for A values of 
0, 20 and 100 meV. Graphs (a) and (b) describe the valence 
band {v\) at F and its backfolded version at L, respectively. 
Graphs (c) and (d) describe the conduction band C3 at F 
(where it follows the top of the original valence band) and L 
respectively. Graphs (e) and (f) describe the conduction band 
ci at F and L, respectively (where it is the original conduction 
band) . 



C3, considering that the Fermi distribution will weaken 
the SW of this band on the theoretical dispersion. At 
T — 65K (Fig. [S] (b), right), the valence band in the 
experimental intensity map shifts to higher binding en- 
ergies, in agreement with the theoretical dispersions cal- 
culated for an order parameter A = 75 meV (Fig. [S] (b), 
left). Moreover, on the experimental map, some inten- 
sity emerges just below Ep (not seen in the false colour 
map, see Ref<^), revealing a dispersive band. Altough 
it does not appear directly in the corresponding calcula- 
tion using 5-peaks, it is reproduced if a finite 30 meV line 
broadening (lifetime) is introduced. In other words, this 
dispersive intensity comes from the combined tails of the 
maximum of the valence band vi (located in the occupied 
states) and of the minima of the backfolded conduction 



9 



Theory | Y | Experiment 



k„(rM) 

-0.2 0.0 0.2 




k„(rM) 

-0.2 0.0 0.2 



T=250K 








f 




1 



T=65K 



it 



0.1 

0.0 

-0.1 7" 

m 

-n 

-0.2 < 

-0.3 

-0.4 
0.1 

0.0 

-0.1 7^ 

m 

-n 
'to 

-0.2 S 

-0.3 

-0.4 



-0.2 



k„(rM) 



0.0 

k„(rM) 



0.2 



FIG. 6: (Color online) Comparisons between theoretical and 
experimental {hv = 31 eV) electronic structures at F. (a) 
The theoretical bands have been calculated for A = 25 meV 
and the experimental intensity maps are taken at T = 250K. 
The continuous black lines highlight the Se4p-derived bands 
not considered in the model, while the dashed white line in- 
dicates the valence band corresponding to vi. (b) The the- 
oretical bands have been calculated for A = 75 meV and 
the experimental intensity maps are taken at T = 65K. The 
dashed black lines indicate the backfolded conduction band 
ci which carries a small non-zero SW (Fig. [5](e)). 



band C3 (located in the unoccupied states). 

Fig. [7] presents comparisons between the theoretical 
and experimental electronic structures at the Brillouin 
zone boundary. The experimental intensity maps at 
T = 250K (Fig. [7] (a)) and T = 65K (Fig. [7](b)) are com- 
pared to calculated bands with A = 25 meV and A = 75 
meV, respectively. At a photon energy of 31 eV, initial 
states between M and L are probed (in a free electron fi- 
nal state picture, we are not probing the BZ exactly at L). 
For simplicity, we will continue using the L notation be- 
low. Substantial changes (explained by excitonic effects) 
in the spectra are evident. At T = 250K (Fig. [7] (a)), 
on the experimental side (right), the conduction band, 
highlighted by the dashed white line, is well described by 
the theoretical band ci (left). It also exhibits a decrease 
of intensity at its bottom, as predicted by our model (see 
Fig. [5] (f)). Moreover it looks like a band having a large 
width and some surprisingly high intensity persists far 
from its centroid, even at binding energies of about 300 
meV. These features can be well explained by the theo- 
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FIG. 7: (Color online) Comparisons between theoretical and 
experimental {hv — 31 eV) electronic structures at the Bril- 
louin zone boundary (see text), (a) The theoretical bands 
have been calculated for A = 25 meV and the experimental 
intensity maps are taken at T = 250K. The dashed white line 
indicates the conduction band corresponding to ci. (b) The 
theoretical bands have been calculated for A = 75 meV and 
the experimental intensity maps are taken at T = 65K. 



retical bands (Fig. [3 left) at A = 25 meV. Due to the 
finite width of real bands, the conduction band ci merges 
with its close (backfolded) neighbours vi or C3 (depending 
on the position along AL), resulting in a band broaden- 
ing. The residual intensity at high binding energies is 
explained with the populated branches of the backfolded 
valence band vi (see Fig. [5](b)). Indeed, looking care- 
fully at the lower part of the experimental map (Fig. [7] 
(a), right), one sees that the residual intensity is larger 
away from L (this is confirmed by momentum distribu- 
tion curves, not shown here). At T = 65K (Fig. [7] (b)), 
on the experimental side (right), a strong signature of the 
CDW appears. The valence band is backfolded at L with 
a high SW. In fact, even a second Se4p-derived valence 
band (indicated by the lowest lying black line in the right 
graph of Fig. [6] (a)) participates to the backfolding at L. 
On the theoretical side, this situation (considering only 
the topmost Se4p-derived valence band) is reproduced 
with an order parameter of A = 75 meV. The backfolded 
valence band vi is well separated from the conduction 
band ci, as in the experiment. It can be clearly seen in 
the corresponding experimental intensity map that the 
conduction band ci does not shift with an increasing or- 
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der parameter (see Fig. [3]) and looses more SW at its 
bottom (see Fig. [5](f)). The backfolded conduction band 
C3 is too far away from Ep in the unoccupied states to be 
measured by ARPES. In the model, at L, the intensity of 
the backfolded valence band vi is lower than that of the 
conduction band ci for high values of the order parame- 
ter, corresponding to a well settled CDW phase (see Fig. 
[5] (b) and (f)). In the ARPES measurements presented 
here, this intensity relation is reversed, as can be seen on 
the right graph of Fig. [7] (b) . The precise reason for this 
matter remains unclear. It can be due to the fact that 
we consider only the topmost valence band in our model. 
Indeed, a second backfolded valence band appears at L 
in the low temperature measurements (Fig. [21(b), right). 
Moreover, we have noticed that this intensity relation be- 
tween the original and backfolded bands can change from 
one sample to another, or even depends on the quality of 
the cleaved surface. Further investigations are needed to 
understand this issue. 



From previous ARPES data^, we extracted an over- 
lap between the valence and conduction bands of 30 
meV. Density functional theory calculation for 1 T'-TiSe2 
within DFT^^ predicts a semimetallic system with an 
overlap of about 800 meV. However this is not completely 
incompatible with our measurements. Indeed, it appears 
that at room temperature excitonic effects are already 
present and the valence band is shifted to higher binding 
energies, below the minimum of the conduction band (as 
we see in Fig. [3] for A ^ meV). It must be emphasized 
that this will be the case, irrelevant of the position of 
the valence band, as long as the valence band maximum 
is above the conduction band minimum in the normal 
phase. In other words, switching on excitonic effects in a 
system having in the normal phase a bandstructure sim- 
ilar to that predicted by DFT would produce room tem- 
perature dispersions very similar to what we measure. 



D. Further discussions 

In the model described in section [III the chemical po- 
tential was not explicitly calculated since it was defined 
as the zero energy of the dispersions. To verify whether 
the chemical potential shifts when the system enters in 
the CDW phase, we have computed the electronic density 
for the renormalized bands, taking into account their dis- 
persion over the whole BZ. Due to the parabolic approx- 
imation of the band dispersions around their extrema, 
we only took into account electrons having an energy up 
to 0.5 eV below Ep. If we keep the chemical potential 
at fj, = eV, a decrease of about 35% of the electronic 
density results between the normal phase and the CDW 
phase with an order parameter A = 75 meV. This dis- 
crepancy is reduced to zero if we shift the chemical po- 
tential by + 60 meV. This result can be understood 
with the SW transfers depicted on Fig. 01 At F, when 
going from the normal to the CDW phase (from Fig. [4] 
(a) to (c)) we loose 13% of the SW of the normal phase 
(integration of the SWs of Fig. [5] (a) and (c)) account- 
ing for the Fermi distribution. At L (from Fig. 2] (a) 
to (c)), after the CDW transition, we loose 40% of SW 
in the conduction band and we acquire 18% of SW in 
the backfolded valence band (integration of the SWs of 
Fig. [S](b), (d) and (f)). Thus, considering only the high 
symmetry directions for illustrative purposes, this results 
in the 35% of SW missing when going from the normal 
to the CDW phase, which can be recovered by slightly 
raising the chemical potential (which affects mostly the 
conduction band ci). 

This shift would be measurable in photoemission. For 
having a detailed knowledge of the chemical potential, 
one needs to perform precise ARPES measurements over 
a wide range of temperatures. However, this is beyond 
the scope of this article and will be studied in the future. 



IV. CONCLUSIONS 

Recently intensity maps calculated within the exci- 
tonic condensate phase model have been compared with 
ARPES data of lT'-TiSe2. Strikingly good agreement 
gave strong evidence for excitonic condensation as the 
driving force of the CDW transition^. In the present 
paper, we have presented the theory of the excitonic in- 
sulator model generalized to the three dimensional case 
of 1 T-TiSe2 with anisotropic band dispersion. From 
the Green's functions of the valence and conduction 
bands, we computed the corresponding spectral func- 
tions needed to generate photoemission intensity maps. 
The mathematical treatment is similar to BCS theory 
and the deduced order parameter in the low temperature 
phase describes the intensity of condensating electron- 
hole pairs (excitons). These pairs are created by the 
electron-hole interaction between the valence band at F 
and the conduction bands at L. As a natural consequence 
of the non-zero momentum of the excitons, this produces 
band backfoldings between F and L which thus tend to 
be equivalent as the order parameter increases. It must 
be emphasized that the CDW produced by this model is 
of purely electronic origin and that the spectral weights 
transferred between the original and backfolded bands 
are large (see reference^ for a more complete discussion 
of this subject). While no real gap opens at the Fermi 
energy, it is notably shown that the valence band (origi- 
nal and backfolded) is shifted in a quasi-linear manner to 
higher binding energies as the order parameter increases. 
Such a behaviour could offer a direct way to extract the 
temperature dependence of the excitonic order parame- 
ter. Finally, the present paper treats only the low tem- 
perature condensation phase of the excitonic condensate. 
However, room temperature measurements indicate that 
strong excitonic fluctuations prevail far above Tc- Their 
theoretical and experimental study promises an interest- 
ing extension of this work. 
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